"Leading blob" model in a stochastic acceleration scenario: the 
case of the 2009 flare of Mkn 501 
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ABSTRACT 

Evidence for very hard, intrinsic 7-ray source spectra, as inferred after correc- 
tion for absorption in the extragalactic background hght (EBL), has interesting 
imphcations for the acceleration and radiation mechanisms acting in blazars. A 
key issue so far has been the dependance of the hardness of the 7-ray spectrum 
on different existing EBL models. The recent Fermi observations of Mkn 501 
now provide additional evidence for the presence of hard intrinsic 7-ray spectra 
independent of EBL uncertainties. Relativistic Maxwellian-type electron energy 
distributions that are formed in stochastic acceleration scenarios offer a plausible 
interpretation for such hard source spectra. Here we show that the combined 
emission from different components with Maxwellian-type distributions could in 
principle also account for more softer and broader, power law-like emission spec- 
tra. We introduce a "leading blob" scenario, applicable to active flaring episodes, 
when one (or few) of these components become distinct over the "background" 
emission, producing hard spectral features and/or hardening of the observed 
spectra. We show that this model can explain the peculiar high-energy charac- 
teristics of Mkn 501 in 2009, with evidence for flaring activity and strong spectral 
hardening at the highest 7-ray energies. 

Subject headings: BL Lacertae objects: general - BL Lacertae objects: individual 
(Mkn 501) - diffuse radiation - galaxies: active - radiation mechanisms: non thermal 
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Introduction 



Blazars are prominent 7-ray emitters with many peculiar features, the origin of 
which are currently discussed intensively in the general context of their multi-wavelength 
properties. During the last years two distinct features have attracted special attention: 
variability detected at very short timescales and the origin of very hard, intrinsic 7-ray 
source spectra when absorption in the Extragalactic Background Light (EEL) is taken into 
account. Although the observed TeV spectra of these sources are steep, their de-absorbed 
(EBL-corrected) spectra appear intrinsically hard. Current uncertainties on the EBL flux 



level and spectrum (cf. 



Primack et al. 



2011 



for a recent review) introduce difficulties 



in defining how hard the absorption-corrected source spectra are. However, in some 
characteristic cases, like for the distant blazars lES 1101-232 {z = 186) and lES 0229-1-200 
{z = 0.139), the emitted spectra still tends to be very hard, with an intrinsic phot on index 



r < 1.5, even when corrected for low EBL fiux levels ( Aharonian et al. 



2006 



200'i . 



The recent Fermi detection of variable 7-ray emission from the nearby {z = 0.034) 
TeV blazar Mkn 501 in 2009 now removes this point of uncertainty, providing strong 
evidence for hard intrinsic 7-ray source spectra independently of que stions related to t he 



20111 ) 



EBL. As already indicated in the original Fermi paper on Mkn 501 (lAbdo et al. 
the spectrum above 10 GeV seems to become much harder during a (~ 30-day) flaring 



Neronov et al. 



fl201lh shows 



state. A recent, independent analysis of the same data by 
that the (10 GeV-200 GeV) flare spectrum could be as hard as F ~ 1.1. While the Fermi 
collaboration did not comment much on the possible origin of the hard flare spectrum. 



Neronov et al. 



(120111 ) put forward the hypothesis that the hard spectrum flare could result 
from an electromagnetic cascade in the intergalactic medium, provided that the strength of 
the intergalactic magnetic fleld is smaller than 10"^^ G and that primary 7-rays with 100 



TeV can escape from the central compact region. 
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The noted small photon indices are generally not easy to achieve in standard leptonic 
scenarios, i.e., synchrotron self-Compton (SSC) or external Compton (EC) models, because 
radiative cooling tends to produce particle energy distributions that are always steeper than 
dN/dE cc irrespective of the initially injected particle spectrum. The corresponding 
TeV photon index would then be F > 1.5. Moreover, as the emission from these objects 
peaks at very high energies where suppression of the cross-section due to Klein-Nishina 
effects becomes important, even steeper intrinsic photon spectra are to be expected. 

Hence, at first glance, hard 7-ray spectra might be somewhat easier to achieve 
in hadronic scenarios like the proton synchrotron model. For a proton spectra of, 
e.g., dNp/dE = AE''^ exp[—{E / Eq)/^] the 7-ray spectrum would be E^dN^/dE^ oc 
j/-('5-i)/2 gxpj— £i//z/g)^/'^+2]. Thus at energies below the synchrotron cutoff hu^ < 3005 GeV 



(lAharonian 



2000) very hard 7-ray spectra may occur. For the magnetic field B ~ 100 G 
required in such models, the protons responsible for 7-rays below huc are not effectively 
cooled. Hence, for s < 2, the 7-ray spectrum could be harder than uF,^ oc z/^/^, with the 
hardest possible z/F^ oc u^^^ in the case of pile- up type proton dis tributions that migh t be 



realized in, e.g., converter-type particle acceleration mechanisms (jPerishev et al. 



20031 ). 



On the other hand, in the context of leptonic models we know that very narrow 
electron distributions are able to produce hard 7-ray source spectra. These can be either 



power-law distributions with large value of the m inimum cut-off, provide c 



field is sufficiently small to avoid radiative los ses (IKatarzvhski et al. 



20091 ) . or provided adiabatic losses dominate (ILefa et al. 



2007 



the magnetic 



Tavecchio et al 



201 ll ). Alternatively, relativistic 



Maxwell-like distributions formed by a sto chastic accelerat ion process that is balanced 



by radiative losses can be a viable option (iLefa et al.ll201ll ). In both cases VHE spectra 

1 /3 

as hard as E^dN/dE^ oc E^ for the SSC case, and E^dN/dE^ oc E^ for EC models 
can be generated. Maxwellian-type particle distributions are especially attractive for the 
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interpretation of the inferred hard 7-ray source spectra because they resemble to some 
extent mono-energetic distributions (the hardest possible injection spectra). An important 
question that arises however then is, whether such distributions can also account for more 
softer and broader photon spectra. 

Interestingly, the hard high-energy component of Mkn 501 emerged in a flaring state 
in May 2009, during which little flaring activity was detected at energies below 10 GeV. 
This suggests that the flare is related to an emission zone that is confined both in space 
(compact) and time (short, on a timescale of ~ 30 da ys). Additionally, there is no evidence 



for a similar, simultaneous increase at X-ray energies (lAbdo et al. 



201l[ ). 



Here we propose a leptonic multi-zone scenario that can accommodate softer emission 
spectra as well as flaring episodes with hard spectral features, like the one observed in 
Mkn 501. In this scenario, the observed radiation comes from several emitting regions 
("blobs"), in which electrons are accelerated to relativistic energies through stochastic 
acceleration, forming pile-up distributions. All blobs are considered to have similar 
parameters apart from the characteristic energies ("temperatures") of their Maxwell- 
like distributions. We note that a "multi-blob" scenario, with power-law components 
characterized by different trajectories (viewing angles ) has al so been proposed i n the 



20081). What 



past to account for the TeV emission from non-aligned AGNs (iLenain et al. 
distinguishes the model introduced here is that, due to acceleration and losses, each 
component is considered to only carry a narrow (pile-up) particle distribution, with broad 
spectra being formed by an ensemble of components. We show, for example, that in the case 
where the total energy of the particles is the same for all blobs, the combined emission leads 
to a spectrum very similar to the one arising from a power-law distribution dN/dE cx 
with index p = 2. On the other hand, distinct spectral feature can appear once a leading 
zone dominates. This could happen, for example, if (a) the value of the temperatures 
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changes from blob to blob, and/or if (b) the energetics of a single blob changes. 



A multi-zone scenario 



Let us consider regions in which electrons are stochastically accelerated (e.g., 
by scattering off randomly moving Alfven waves) up to energies where acceleration is 
balanced by synchrotron or inverse Compton (Thomson) losses. Their steady- state energy 
distributions n('y) oc l^ fil) then take on a relativistic Maxwellian-type form (jSchlickeiser 



1985 



Aharonian et al. 



198' 



0; i = 1, .., A^) with exponential cut-off Lorentz factor 



7c 



bD, 



1/b 



(2) 



and normalization factor Ai. Here, Ps refers to energy losses due to synchrotron radiation 
given by 

fin _ ZL _ _ _ 

(3) 

DnV^ -\ with 



dt 3 

and the constant Dq is given by the diffusion coefficient Dp = ^{^^ 



Va 



B 



/Attp 



the Alfven speed and t = X/c cc p ^, b > 1, the mean scattering time ( iLefa et al. 



20111 ). The shape of the exponential cut-off (characterized by the parameter b) is related to 
the turbulence wave spectrum W{k) oc fc""^ as b = q. Note that if the particle distributions 
would be shaped by inverse Compt on cooling in the Klein-Nis hina regime, a smoother 



exponential cut-off is expected (e.g.. 



Stawarz fc Petrosian 



20081 ). 



For simplicity, we consider below the situation where all blobs have similar properties 
(e.g., magnetic field strength, linear size, Doppler factor) but different values for the 
characteristic energy 7c^. In the case of synchrotron losses and scattering off Alfven waves, 
7c. depends on the magnetic field and the bulk density of the flow. Thus, the modification 
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in 7c. that we assume for each blob might be related to a non-homogeneous bulk flow. 
The total energy density Ei that the relativistic particles gain through scattering off Alfven 
waves is calculated to be 

Ei= 'jni{'y)d'ymec'^ = AinieC^ -^T[-, — ] , (4) 
Ji 7^. 

where r[a,2;] is the incomplete F function. We can now express the normalization factor as 
a function of the "temperature" (7c J and the energy Ei, 

which decreases as the cut-off energy increases, with a dependency well approximated by 
Ai oc 7~^ for all values of the coefficient h of interest. This is easily seen if we change the 
lower limit of the previous integration to 0, in which case the result becomes 

A- = (Q) 
' m,c'^lT[{A + b)/b]' ^' 

with r[z] denoting the F function. 

The combination of the above electron distributions can lead to power-law-like particle 
distributions if the temperatures of the different components do not differ significantly and 
if the components contribute equally to the overall spectra. The power-law index then 
essentially depends on the amount of energy given to the non-thermal particles in each 
blob, i.e., on how the total energy E^ of each component scales with temperature (7c J • 
Steep spectra may arise if, e.g., the low temperature components dominate whereas harder 
spectra may occur if more energy is contained in the high-temperature blobs. 

In Fig. [T] an example for the total differential electron number is shown assuming 
Ei = constant. For the plot, an exponential cut-off index 6 = 2 and = 4 have been 
chosen. The temperatures are equally spaced on logarithmic scale. Then, the total energy 
distribution approximately forms a power-law dNf./d'-f oc 7"'* with index s ~ 2 between 
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the minimum and maximum temperatures. As discussed above, this "special" value of the 
power-law index results from the assumption that all blobs have the same total energy, 
so that 7^(iAi"e/(i7 is the same for each zone. This can be demonstrated more formally by 
looking for the "envelope", i.e., the mathematical function that describes the curve which 
is tangent to each of the curves rii in the {dNe/d'-f,'-fc) plane at some point. This function 
approximates the sum of the energy distribution of the different components as >> 1 
and gives the characteristic behavior of the total distribution above the minimum and below 
the maximum temperature. It can be found by solving the set of equations 

F(7,7,) = 0, ^,F{j,j,) = 0, (7) 

where F(7,7c) = dN^./d'y — v4(7c) 7^ exp[— (7/7C)'']. If we assume that the energy in 
non-thermal particles is the same for each component, then we find 

dNe/d'J = C'{b) 7"^ 7c,min < 7 < 7c,max , (8) 

where c'{b) = 4:{4:/b)^/^e-^/^/meC^T[^] (cf. Fig. [1]). On the other hand, if particle 
acceleration to higher energies goes along with a decrease in total energy, e.g., E cx l/7c, 
then steeper power law spectra can appear, i.e., oc 7"^. Note that in all these cases, 
radiative cooling is already taken into account. 

Note that changing the total energy with in each component could also be 
interpreted as changing the number of contributing blobs as a function of 'Ja, assuming that 
each component has the same energy. This could be formally accommodated by introducing 
a statistical weight Wi, so that the overall spectrum is expressed as 

i=N 

dNjd^ = Y,wM^). (9) 

i=l 

Harder spectra may then occur, for example, if more blobs with higher temperatures exist 
and vice versa. Hence, a conclusion similar to the above can be drawn, once the statistical 
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weights vary correspondingly with temperature (7cJ- In a continuous analogue, we may 
write 

dNjd-i = V wM-i) ^ / (T)n(7, T)dT (10) 

where W{T) is the spectrum of the number of components per temperature. Since 
Maxwellian-type electron distributions behave, to some extent, like mono-energetic ones, 
the total energy distribution between the minimum and the maximum temperature mimics 
the spectrum of the number of the blobs W{T) (provided W{T) does not rise quicker than 
T^). Thus, in principle a variety of spectra may arise, depending on the choice of W{T). 
Conversely, observations of extended power law-like energy distributions then impose 
constraints on how W{T) of a source can vary with temperature. 

The SSC spectrum, arising as the sum of the different Maxwell-like distributions 
of Fig. [H is shown in Fig. |2l The synchrotron flux resembles the flux that would be 
emitted by a power-law particle distribution of index 2. Between the frequencies related 
to the minimum and maximum temperatures, v^i^ cx -B7cmin ^'^d t'max -^Tcmax 
exhibits a F,^ oc z/^^^^ behavior. For v < z/mm it follows the characteristic F^, oc u^^^ 
synchrotron emissivity function, w hile for > Uma^ the exponential cut-off b ecomes 



smoother oc exp [— (z^/z^max) ^+''] ( Fritz 



1989 



Zirakashvili fc Aharonian 



20071 ). 



In the present model, the electrons in each blob are considered to only up-scatter their 
own synchrotron photons and not the ones emitted from the other blobs. The photon fields 
produced by the other components therefore do not contribute to the emitted Compton 
spectrum of each blob. In the Thomson regime then the up-scattered photon spectrum 
again approaches a power-law behavior similar to the synchrotron one, i.e. F^, oc z/^^/^. 
Once Klein-Nishina effects become important, the suppression of the cross-section makes 
the high-energy spectrum steeper, as expected. In the case of discrete zones, where the 
particle distributions and synchrotron photons are almost mono-energetic, this happens 
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at (intrinsic) energies greater than ^ci{B /Bcr)^1. > 1, where Bcr — mlc^/{eh). Here, 
{B/ Bcr)l1- is the peak energy of the synchrotron photons emitted by the i-th blob with 
temperature 7c.. Below u oc 1c,minBlc,min Compton flux reveals the characteristic 1/3 
slope, reflecting the low-energy synchrotron spectrum. In the Klein-Nishina regime, the 
exponential VHE cut-off mimics the shape of the electron cut-off, Fi, oc exp [— (i^/^'max)''], 
and is steeper compared to the synchrotron spectrum. 

3. The origin of hard 7-ray flare spectra 

Once a single component becomes dominant in the overall emission, as naturally 
expected for a flaring stage, hard spectral features can arise. This is more evident in the 
Compton part of the spectrum as (in the Thomson regime) the separation of the VHE 
peaks scales as ~ 7^ and is greater than in the synchrotron case (~ 7^). The energetics 
of such a leading component, which is responsible for an observed flare, could change for 
different reasons: The total (intrinsic) energy offered to the accelerated particles or/and 
the temperature of the distribution could increase, for example, due to changes in the bulk 
flow properties or due to an increased injection of seed particles that undergo stochastic 
acceleration. Another possibility concerns an increase of the Doppler factor. Already a 
shght change of the hne-of-sight angle to the observer during the propagation of a blob, 
could lead to the observation of a hard flare without an accompanying change of the 
intrinsic energetics of the components. 

The aforementioned possibilities can be applied to explain the high-energy flare of 
Mkn 501 observed in 2009. To illustrate this, we assume for the "low state" that the total 
energy of each component drops as the temperature increases in order to account for a 
steeper than 1.5 spectrum. The scale assumed is £^ oc ^c^^^. For the flaring state the 
normalization of the two components with the highest temperatures is increased by a factor 
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of ~ 2. Their temperature is also slightly increased. This leads to a hard high-energy flare 
above 10 GeV. Below 10 GeV the emission stays almost constant. Other parameters are 
kept the same for all blobs (i.e., magnetic field B = 0.1 G, blob radius R = 10^^ cm and 
Doppler factor 6 = 30). For this set of parameters, the synchrotron flux is lower than the 
one observed, implying that X-ray emission would have to come from a different part of the 
jet. This is consistent with the fact that little flux variation has been observed in X-rays. 



4. Conclusion 

Narrow energetic electron distributions, like relativistic Maxwellian-type ones, can 
successfully explain the very hard intr insic 7-ray spect ra that arise in some sources once 



EBL absorption is taken into account (ILefa et al. 



201l|). 



Here we have demonstrated that the superposition of emission from such distributions 
could also accommodate more softer and broader 7-ray spectra. To show this, a multi-zone 
scenario was considered in which particles are accelerated through a stochastic acceleration 
process balanced by radiative (synchrotron/Thomson) losses in multiple zones characterized 
by different temperatures (i.e., achievable maximum electron energies). For the parameters 
examined here, particle escape can be neglected, and the particle distribution in each zone 
essentially takes a Maxwellian-type shape. Under reasonable conditions, the resultant 
overall (combined) particle energy spectra then approaches a power-law particle distribution 
dNf./d'-) oc 7~'^ over the energy range corresponding to the lowest and the highest 
temperature, with power index s only depending on how the total energy (in non-thermal 
particles) in each component scales with temperature (cut-off Lorentz factor 7c). In the 
case where all parameters, apart from the temperature, are kept constant (in particular 
the total energy in each component), the resultant power index approximates s ^ 2. For 
similar magnetic fields and Doppler factors, softer/harder 7-ray spectra could arise when 
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the lower/higher temperature components dominate. 



In this scheme, the dominance of one (or a few) of the radiating components could 
lead to a flaring state during which hard spectral feature become apparent. This leading 
component might increase its luminosity for different reasons, e.g., due to a change of the 
Doppler factor or the injected energy. As shown above, such a scenario can account for the 



2009 flare in Mkn 501 d uring which a strong hardening o : 



GeV has been observed (lAbdo et al. 



2011 



Neronov et al 



the emission spectra above 100 



2011 



Mkn 501 is indeed known 



to be a source where detailed temporal and spectral modeling has provided evidence for the 
contributions of diff erent components (such a s a steady X-ray component plus a variable 



SSC component, see 



Krawczynski et al. 



20021 ). 



While in the present work an SSC approach has been employed, similar features are to 
be expected in external Compton scenarios. In the latter case, an even stronger spectral 
hardening up to Fy (x u may occur, while in the SSC case this is limited to Fy cx z/^/^ 



( iLefa et al 



20111 ) 



In summary, the "leading blob" model introduced here can explain in a natural way 
both the quiescent and flaring state in Mkn 501. It will be interesting to check to what 
extent this also applies to other sources. 
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2 2.5 3 3.5 4 4.5 5 5.5 6 

Logy 



Fig. 1. — The total electron energy distribution from = 4 blobs with different temperatures 
(7c) but the same total energy Ei tends to establish a power-law dNe/d'^ oc 7"^ between the 
minimum and maximum characteristic energies 7c. Below 7c,miru dNf.jd'y oc 7^, whereas for 
7 > 7c,max one has dNf.jd^ oc e"'^'''/''''^'™''''^^. Here, the exponential cut-off index is 6 = 2 and 
other parameters used are i? = 0.1 G and Ei = 2 x 10^^ erg. 
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Fig. 2. — Resultant SSC emission from the combination of different components with 
Maxwelhan-hke electron distributions. The synchrotron flux exhibits a power-law behav- 
ior, oc z/~^/^, approximately between the energies related to the minimum and maximum 
temperature. The same holds for the Compton flux in the Thomson regime, while in the 
Klein-Nishina regime the spectrum becomes steeper. Doppler factor 6 = 30 and cut-off index 
6 = 2 have been used. Other parameters are R = 3 x 10^^ cm, B = 0.1 G, Ei = 2 x 10^^ erg. 
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Fig. 3. — Resultant SSC emission from the sum of different (A^ = 4) blobs for the "low" 
(dashed-line) and "flaring state" (blue line) of Mkn 501. The total energy given to the par- 
ticles scales as Ei oc 'y^^^^ and for the low state the temperatures are log(7c) = 3.9,4.14,4.6 
and 5.2, respectively. Parameters kept constant are magnetic field i? = 0.1 G, blob radius 
R = 10^^ cm, cut-off index 6 = 3 and Doppler factor 6 = 30. For the flaring state, the two 
blobs with highest temperatures are assumed to be enhanced by a factor of ~ 2 with their 



temperatures slightly increased (to 4.74 and 5.3, respectively). Belo w 10 GeV, t 



almost constant wit h respect to the low state. For data points, see 



Neronov et al. 



2011 



Abdo et al. 



le flu x is 
2011 and 
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Fig. 4. — Same as figure [3] but assuming the flaring state to occur due to a change of the 
Doppler factor of the two components with the highest temperatures from 5 = 30 to 5 = 40. 
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Log E [eV] 

Fig. 5. — Same as figure H] with the synchrotron part of the spectrum included. The X-ray 
regime is considered to be dominated by emission from a different part of the jet. Thus, 
during the flare no/little variability at X-ray energies would be observed as long as the 
synchrotron contribution from the "flaring" components does not exceed the measured X- 
ray data. 



